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An optimization procedure determining tiie ideal configuration at the microstructural level of 
ferroelectric (FE) materials is applied to maximize piezoelectricity. Piezoelectricity in ceramic FEs 
differ significantly from that of single crystals because of the presence of crystallites (grains) pos- 
sessing crystallographic axes aligned imperfectly. The piezoelectric properties of a poly crystalline 
(ceramic) FE is inextricably related to the grain orientation distribution (texture). The set of combi- 
nation of variables, known as solution space, which dictates the texture of a ceramic is unlimited and 
hence the choice of the optimal solution which maximizes the piezoelectricity is complicated. Thus 
a stochastic global optimization combined with homogenization is employed for the identification 
of the optimal granular configuration of the FE ceramic microstructure with optimum piezoelectric 
properties. The macroscopic equilibrium piezoelectric properties of polycrystalline FE is calculated 
using mathematical homogenization at each iteration step. The configuration of grains characterised 
by its orientations at each iteration is generated using a randomly selected set of orientation dis- 
tribution parameters. The optimization procedure applied to the single crystalline phase compares 
well with the experimental data. Apparent enhancement of piezoelectric coefficient dss is observed 
in an optimally oriented BaTiOa single crystal. Based on the good agreement of results with the 
published data in single crystals, we proceed to apply the methodology in polycrystals. A config- 
uration of crystallites, simultaneously constraining the orientation distribution of the c-axis (polar 
axis) while incorporating ab-plane randomness, which would multiply the overall piezoelectricity in 
ceramic BaTiOa is also identified. The orientation distribution of the c-axes is found to be a narrow 
Gaussian distribution centred around 45°. The piezoelectric coefficient in such a ceramic is found 
to be nearly three times as that of the single crystal. Our optimization model provide designs for 
materials with enhanced piezoelectric performance, which would stimulate further studies involving 
materials possessing higher spontaneous polarization. 

PACS numbers: 77.65.Bn,77.65.-j,77.80.-e,02.60.Pn 



I. INTRODUCTION 

Ferroelectrics (FEs) have plethora of applications from 
actuators to sensors and from ultrasonic generators to 
energy-harvesting devices due to its piezoelectricity (or 
electric-field-induced strain)ir— . The piezoelectric coeffi- 
cients d, which quantify the piezoelectric strain against 
an applied electric field, are the ubiquitous figures of 
merit in such applications^i^i^. Optimization of these pa- 
rameters in single crystals and nanoparticles compounds 
with the existing deleterious effects such as twinning 
and depoling though many of them exhibit enhancement 
of piezoelectricity in certain nonpolar directions^i^ii. 
(Here nonpolar in the sense that a direction other than 
the spontaneous polarization direction of the crystal). 
Nonetheless, FEs in the polycrystalline form are pre- 
ferred over the single crystals in engineering owing to the 
ease in manufacturing and the compositional modifica- 
tions of polycrystals®. NoncoUinear polarization rotation 
has been proposed as the possible origin of high piezo- 
electric response in FE single crystals^i^"— . Furthermore, 
phenomenological thermodynamic theory addressing the 
piezoelectric anisotropy relates it to the flattening of the 
free energy function in certain nonpolar direction s 

As-grown polycrystalline FE is an aggregate of single 
crystalline grains with randomly oriented (spontaneous) 
polarizations^^. The spatial configuration of crystallo- 



graphic grains and their orientation distribution (tex- 
ture) impact the piezoelectricity exhibited by conven- 
tional as well as new generation FEsi. The randomness 
in polarization-vector orientation renders the resultant 
piezoelectricity of the material to be marginal or zero. 
Although the resultant polarization is zero for as-grown 
polycrystal, an overall piezoelectricity can be enabled by 
the application of an external electric field, called poling 
field, though all the grains can never align perfectly^'*. 

The aggregate texture of an unpoled polycrystal would 
have a uniform random distribution of orientatio n^^' 
With the strength of the poling field increases, we as- 
sume the nature of the grain orientation distribution to 
become Gaussian (normal) "'^^d'' . Piezoforce microscopic 
image studies of Pb(Zr,Ti)03 (or PZT) films leads also to 
a quantitatively similar kind of distribution of domains^-. 
Spherical harmonics based orientation distribution func- 
tion (ODE) is also used to model the texture of poly- 
crystalline ferroelectrics—-"'^ . The texture analysis of 
PZT by synchrotron X-ray diffraction reveals that the 
diffraction peak intensity ratio i?{2oo} shows a strong de- 
pendance on sin^ 9, where 9 represents the orientation 
angle between the plane normal and the polar axis of the 
materiaP^. Moreover, the domain distribution has been 
treated statistically to be in an orientation distribution 
characterised with March-Dollase function in X-ray and 
neutron diffraction studies of piezoelectric materials^^. 
fn this work, Rogan et ali^ has done in situ neutron 
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diffraction study of tlie polycrystalline PZT wliile me- 
clianical loading. Tliis work provides insiglits into tlie fer- 
roelastic switching of domains in FEs and uses a March- 
DoUase function to model the evolved texture during 
loading. 

Electrical and mechanical loading (or poling) have 
limitations in realising smooth saturation textures due 
to the restrictions imposed by the crystallographic 
symmetry2ii2^. Many domains cannot be reoriented due 
to a complex set of internal stresses and electric fields 
in grains and because some domains will switch back 
after the poling field is removedS.. Yet, the Gaussian 
distribution is found to realize the piezoelectricity espe- 
cially at the boundaries bounded by a random distribu- 
tion (where the piezoelectric coefficients, djiy — ^ 0) and 
a perfect single crystal texture (where dj^, tends to the 
single crystal values) correctly^^. This can overcome, for 
instance, the inadequacies of other distribution functions 
in approximating the case of strong preferred orienta- 
tion, which eventually refiect in the correct evaluation of 
piezoelectricity^!. For instance, March-Dollase function 
models the distribution very efficiently at low and moder- 
ate degrees of non-randomness. But the strong preferred 
orientation cannot be adequately approximated by this 
functional. In the present model we present configura- 
tion of crystallites with orientation distribution fits into 
a Gaussian. Despite this, we are not deterministic about 
the methods (e.g., hot forging, reorientation by poling, 
tape casting) to realise this grain configuration (distri- 
bution). Nor are we attempting to model the evolution 
texture of piezoelectric ceramics in this paper albeit pre- 
dicting the configuration of grains that would optimize 
the piezoelectricity. 

The spatial and orientational randomness of grains can 
be judiciously employed in the design of EE polycrys- 
tals (ceramics) with tailor-made configurations^^. Re- 
cently it is shown that overall piezoelectricity would be 
enhanced by the introduction of either [110]- or [111]- 
oriented grains into a random BaTiOa polycrysta l^^'^" . 
The piezoelectric coefficient 6^33 scales up until a certain 
grain size and further decrease from this point would 
decrease the ^33 in these experiments. Though the ce- 
ramics possess randomly oriented grains, fine domains 
must be introduced to achieve this narrow window of 
high piezoelectric activity. The solution of the optimal 
grain configuration of the polycrystalline EE is hardly 



resolved although much research activity is underway in 
ferroelectrics. This is mainly due to the vast number 
of possible configurations available at hand albeit the 
EE ceramics are easy to manufacture. In this work we 
would identify an optimum configuration of grains in the 
microstructure of a polycrystalline EE material with an 
enhanced piezoelectricity from both a textured ceramic 
and a single crystal. Moreover, we would discuss on the 
mechanism of polarization rotation in EEs in view of their 
piezoelectric response in the context of the optimization 
results. 

Eor many EE single crystals, the polarization P rotates 
with an electric field applied along a nonpolar direction 
and the polarization rotation in such cases does not nec- 
essarily evolve through a monoclinic phased. The polar- 
ization Pi = di^Tf^, where T is the stress and d piezo- 
electric strain coefficient. Here the Latin indices range 
from 1 to 3 and Greek indices from 1 to 6, owing to the 
Voigt contraction of tensorial indices^. The typical value 
of the piezoelectric coefficient d^s measured for tetrago- 
nal (T) BaTi03 is w 90 pC/N when [001]-poled while 
1^33 w 203 pC/N for [lll]-poled engineered-domain sin- 
gle crystal^i^. The piezoelectric coefficient dsi typically 
falls between -33.4 and -62.0 pC/N for [001]- and [111]- 
poled single domain BaTiOs crystals^^jJ^,. Since tem- 
plated grain growth^! enables the fabrication of textured 
ceramics with a fraction of oriented material, it is impor- 
tant to find out the degree of orientation as well as the 
fraction of aligned grains in a matrix of random polycrys- 
talline EE. Pb-free (Ko.5Nao.5)(Nbo.97Sbo.o3)Oo.3 poly- 
crystal grown using this method is recently shown to pos- 
sess excellent electromechanical properties if it is (001)- 
oriented with a narrow orientation distribution^. Such 
a possibility to fabricate useful polycrystals is one of the 
motivation to a continuum model study, as is shown in 
this paper. 

The orientation of a piezoelectric crystal is modeled by 
a set of Euler angles (0, 9, ip). Euler angles are defined in 
the following ways: first, the crystal is rotated by angle cj) 
around the z-axis, then rotate an angle 6 around the new 
x-axis, and finally by an angle ip around the new z-axis. 
All the rotations are in the counterclockwise direction"^—. 
The matrix of transformation from the crystallographic 
coordinate system to a local coordinate system y is given 
by^ 



cos Ip COS 4> — COS 9 sin (j> sin ?/; cos sin (/) + cos 9 cos (p sin ip sin -0 sin ( 
— sin Ip cos 4> — cos 9 sin (p cos ip — sin ip sin (p + cos 9 cos (p cos ip cos ip sin ( 
sin 9 sin 6 — sin 9 cos 6 cos [ 



(1) 
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The crystal orientation and thus the piezoelectric proper- angles. Even in the case of single crystals it would be te- 
ties would therefore inextricably depends on three Euler dious to survey all the crystallographic orientations and 
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map the corresponding piezoelectric properties. This sit- 
uation demands an optimization procedure to reahze a 
particular objective function in a piezoelectric material. 
[There would be 73'^ {— 389,017) different combinations 
of {(j>, 6, ip) if one chose angles at an interval of 5° between 
±7r in FE single crystals alone]. Nonetheless, the case of 
polycrystals is too complicated as it is constituted ide- 
ally by thousands of single crystalline grains. To arrive 
at an optimum texture of the ferroelectric polycrystal at 
which the material exhibits maximum piezoelectric per- 
formance, a global optimization method has to be em- 
ployed here as well. This is because the piezoelectricity 
depends on the parameters which controls the orienta- 
tion distribution of the grains. Nevertheless, the choice 
of the optimal set of parameters is complicated and it is 
impossible to analyze all possible combinations of the dis- 
tribution parameters or the angles themselves. We chose 
a modified stochastic global optimization technique in- 
corporating a generalised Monte Carlo scheme for this 
purpose. A modified simulated annealing (SA) is quite 
suitable in this respect as the objective function is not 
sensitive to the starting point of the iterative process (the 
so called connectivity, where any state of the system can 
be reached starting from any other state)2i. 

We propose a design strategy based on a continuum 
mechanics and modeling to attain electromechanic fig- 
ures of merit in FEs combined with a modified simulated 
annealing for optimization, in this paper. The model is 
developed for crystals of all classes. In this paper, we 
apply this methodology to optimize the piezoelectricity 
in the classic perovskite, the tetragonal PAmm BaTiOa. 
The piezoelectric performance is quantified through the 
evaluation of the effective electromechanical properties 
of the FE single crystals and polycrystals. We have used 
the mathematical homogenization method^^i^S which ef- 
ficiently characterizes the equilibrium macroscopic elec- 
tromechanical properties of a polycrystalline ferroelectric 
material. 



II. MODEL 



tion of different phases in characterizing both the macro- 
and the micro-mechanical behaviors. (For instance, in 
a polycrystalline material the crystallite manifests as a 
phase in this sense). In homogenization theory the ma- 
terial is locally formed by the spatial repetition of very 
small microstructures (unit-cells), when compared with 
the overall macroscopic dimensions. Further, the mate- 
rial properties are periodic functions of the microscopic 
variable, where the period is very small compared with 
the macroscopic variable. This enables the computation 
of equivalent material properties by a limiting process 
wherein the microscopic cell size is approaching zero. 

The key idea in micromechanical modeling is to 
relate the effective properties of a material to the 
properties of its constituent parts, which may be the 
phases of a composite or the grains of a polycrystal in 
ferroelectrics^Si^^ii^Tii^. A review on the various meth- 
ods for obtaining effective properties of ferroelectrics is 
given in an earlier paper—. The asymptotic analysis and 
homogenization of the piezoelectric medium has resulted 
in the equilibrium macroscopic piezoelectric properties in 
tensor notation as follows; 
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A. Homogenization of ferroelectrics 

The electrical and elastic boundary conditions and the 
orientation of the polarization and permittivity axes of 
grains (in polycrystals) and crystals complicate the mi- 
croscopic analysis of the FEs^i^. Considering these, 
the mathematical modelling (if we describe the rapidly 
varying material properties with equally rapidly varying 
functions) and the numerical analysis of these materials 
will become difficult and sometimes even intractable. To 
put it simply, homogenization of partial differential equa- 
tions (PDEs) has as its main purpose to approximate 
PDEs that have rapidly varying coefficients with equiv- 
alent homogenized PDEs that (for instance) more easily 
lend themselves to numerical treatment in a computer. 
The homogenization method accommodates the interac- 



Here the e, k and C are the electromechanical coeffi- 
cients, viz., piezoelectric, dielectric and elastic stiffness 
coefficients respectively. Also, xVi^^y) is the charac- 
teristic displacement, R{x,y) is the characteristic elec- 
tric potential, ^f{x,y) and iprs{x,y) are characteristic 
coupled functions of the FE unit-cell of size Y, satisfy- 
ing a set of microscopic equations^^i^. 6 is the Kro- 
necker delta symbol. (We have assumed Einstein sum- 
mation convention that repeated indices are implicitly 
summed over throughout this paper). Symmetry requires 
that Bktj = ekji, Kfj = K.'ji and C^^ = C^^. In the 
above Eqs. ©-(g]) ekij{x,y), Kfj{x,y) and C^^[x,y) are 
the electromechanical properties of the single crystallite 
whose collection constitutes the microscopic unit-cell and 
can be described in microscopic coordinates y as 
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G^ipdjqfXkr^pqr 



(5) 



where are the Euler transformation tensors from crys- 
tallographic coordinate system to the local microscopic 
coordinates y-^ . Here the primed moduli are the ones ex- 
pressed in crystallographic coordinate system. Also the 
piezoelectric response is determined along an arbitrary 
crystallographic direction determined by the Euler an- 
gle ((/), 0, ij}) with respect to the reference frame of the 
microstructure, i.e., y. [Here the microstructure refers 
to the representative volume element (RVE) of the fer- 
roelectric polycrystal used for the homogenization] . The 
(resultant) polarizations of the crystallites in an as-grown 
polycrystal are randomly oriented in the lattice space and 
hence require three angles to describe its orientation with 
reference to a fixed coordinate system. Euler angles 0, 
6 and tp can completely specify the orientation of the 
crystallographic coordinate system embedded in crystal- 
lites and thereby the orientation of relative to a fixed 
Cartesian coordinate system. (The superscript H would 
be dropped from the homogenized piezoelectric proper- 
ties given in Eqs. (HI -(HI for brevity from the rest of the 
discussion) . 

After asymptotic analysis a set of microscopic system 
of equations characterising xV{^Ty)^ ^(^J,?/), ^^{x^y) 
and iprs{x,y) is obtained and is solved computation- 
ally. The three-dimensional (3D) numerical model de- 
veloped is implemented in finite element method (FEM) . 
The polycrystal sampling is performed using a unit-cell 
of 14 X 14 X 14 mesh with 21,952 Gaussian integration 
points. The convergence of piezoelectric properties with 
unit-cell size allows us to determine the simulation-space- 
independent, macroscopic piezoelectric properties at var- 
ious distribution of grains. Also the convergence studies 
are necessary as the orientations (Euler angles) of the 
crystallites are chosen randomly from a Gaussian distri- 
bution of angles. A normal random generator delivers 
different sets of numbers each time it is invoked and this 
will affect the consistency of the results. Nonetheless, at a 
14'^ mesh we found that the statistical fluctuations on di^ 
are less than 4% (Ref.[2^. After using the usual approxi- 
mations of FEM, the set of linear equations for each load 
case is obtained where each global stiffness, piezoelectric 
and dielectric matrix is the assembly of each element's 
individual matrix, and the global force and charge vec- 
tors are the assembly of individual force vectors for all 
the elements. 

Full integration (2-point Gaussian integration rule in 
each direction) is used for the evaluation of the stiffness, 
piezoelectric and dielectric matrices and for the homoge- 
nization. As the representative microstructure (unit-cell) 
is expected to capture the response of the entire piezo- 
electric system, particular care is taken to ensure that 
the deformation across the boundaries of the cell is com- 
patible with the deformation of adjacent cells. Thus all 



the load cases are solved by enforcing periodic bound- 
ary conditions in the unit-cell for the displacements and 
electrical potentials. Though the model is general and 
is designed to accommodate all crystalline symmetries, 
we apply it to the case of BaTiOa. The numerical ho- 
mogenization of ceramic BaTiOa is carried out using the 
single crystal data taken from Ref. |32j using the present 
homogenization model. 



B. Optimization of piezoelectricity 

The texture of a polycrystalline ceramic can be quan- 
tified through the distribution parameters corresponding 
to the angle of orientation in space. Euler angles (0, 9, -0) 
can completely specify the orientation of the crystallo- 
graphic coordinate system embedded in crystallites and 
thereby their orientations with respect to a local coordi- 
nate system j/i. Since the aggregate texture for polycrys- 
talline FEs follows a Gaussian distribution, the probabil- 
ity distribution function (pdf) is defined by. 
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about the direction of the electric field, fi and a are the 
parameters of the distribution viz., the mean and the 
standard deviation respectively, a stands for the Euler 
angles (0, 9, ip). Since /i and a decide the scatter of orien- 
tations (of the grains) and for that matter be critical to 
the piezoelectric response of an FE ceramic, they would 
assume the role of design variables of the optimization 
problem. Thus we are aiming to find an optimum set 
of these parameters from a solution space controlled by 
the laws of coordinate transformations from a crystallo- 
graphic coordinate system embedded in the grains to a 
local coordinate system which coincides with the global 
frame of reference. Also, the solution space is bounded 
by distribution parameters ranging from those of uniform 
(in the case of random polycrystal) to those of normal 
distribution (in the case of poled polycrystal). A fairly 
uniform distribution can be achieved by putting standard 
deviation (cr) equals 5 and for a poled ceramic ferroelec- 
tric the a is set near zero. 

Here the objective is to search possible ways of en- 
hancing the piezoelectricity in FE polycrystals. When 
an electric field is applied in a piezoelectric crystal the 
shape of the crystal changes slightly. This is known as 
the converse piezoelectricity. There exists a linear re- 
lationship between the components of the electric field 
vector Ei and the components of the strain tensor Sij 
which describe the change in shape"^^. The piezoelectric 
strain constant d^^, (which is written in tensor form as 
£jk = dijkEi) thus quantifies the piezoelectricity of a FE 
material. In a ceramic FE each of the three Euler an- 
gles 0, 9 and ijj if observed individually falls in respective 
normal distributions thanks to the misalignment of the 
crystallographic axes of the constituent grains. Hence 
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there would be six (3 pairs) parameters altogether viz., 
(/^^jf^); (MSjCTfl) and (11^,0^) quantifying the scatter of 
angles from j/j. The grain distribution parameters cho- 
sen by the optimization algorithm will prompt a normal 
random generator and thereby create a set of Eulcr an- 
gles (0, 6, ip) for each of the grains. These Euler angles 
will dictate the coordinate transformation [as given by 
Eq. ([5])] in the electromechanical property tensors ap- 
pearing in the homogenization equations Eqs. ([2]) - Q. 

Each iteration of the optimization algorithm calls the 
objective function, the effective piezoelectric coefficient 
dj^ . The numerical solution of the coupled piezoelectric 
problem sought using the FEM would be substituted in 
the homogenization expressions. Nonetheless, the elec- 
tromechanical tensor undergo coordinate transformation 
dictated by the Euler angles as depicted in the expres- 
sions in Eqs. ([5]) before it enter the homogenization algo- 
rithm. The FEM used for this study correlates each ran- 
domly oriented grain in a polycrystalline material with 
each element of the finite element mesh. The macroscopic 
piezoelectric coefficients e^^ obtained in the homogeniza- 
tion equation Eq. ^ is calculated first and subsequently 
we calculate = Tlifi=i^i'j^^v^ tiy rnaking use of the 
elastic compliance derived from inverting homoge- 
nized C^^ from Eq. gj. 

The key idea behind simulated annealing is based 
in the Monte Carlo step proposed in the Metropolis 
algorithm^ for simulating the behavior of an ensemble 
of atoms that are cooled slowly from their melted state 
to the minimum energy ground state. The ground state 
or minimum energy state corresponds to the global opti- 
mum we are seeking in material optimization. In order 
to apply this algorithm to a piezoelectric material, we 
must first introduce the notion of system energy. In the 
present setup the piezoelectric coefficient ^33 acts as sys- 
tem energy and we seek its maximization. In order to be 
consistent with our definition of design variables, let 

E{R,) = dssia)^ 

E{Ri) = d^s^a^, fi^,ae, fJ-e,cr^, ^ji,) (7) 

be the surrogate of energy E of a particular configuration 
Ri. Here each set of distribution parameters are selected 
randomly from a € [0,5] and fi S [0,7r/2]. The main 
goal of optimization is to find the ground state(s), i.e., 
the minimum energy configuration(s), with a relatively 
small amount of computation. Minimum energy states 
are those that have a high likelihood of existence at low 
temperature. The likelihood that a configuration, Ri, 
is allowed to exist is equal to the Boltzmann probability 
factor, P{Ri) = exp[—^^^^\, where /cb is the Boltzmann 
constant and T is the temperature, ks is often treated 
as unity for computational convenience**^. 

A control parameter similar to the temperature in 
physical annealing is introduced in optimization which 
will dictate the number of states to be accessed in going 
through the successive steps of the optimization algo- 
rithm before being settled in the minimum energy state 
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FIG. 1: (a) Objective function (piezoelectric coefRcient d^z) 
as a function of temperature in FE single crystal BaTiOa . Tiie 
filled diamond indicates tiie experimental value {= 203 pC/N) 
obtained in Ref. 's^for [lll]-oriented engineered-domain sin- 
gle crystal BaTiOs. The dotted line shows the evolution of dzs 
with iteration. (Here the temperature is a decreasing control 
parameter of the optimization). 




FIG. 2: The piezoelectric coefRcient dss (spheres) derived in 
the optimization as a function of Euler angles; (a) 6 and cj> 
and (b) 6 and 



(the optimum configuration). In single crystals, the de- 
sign variables are the three Euler angles (0, 6, ip) itself. 
The design domain is constrained by — tt < ((j), 9, -0) < tt. 
This condition would be enough to scan the entire crys- 
tallographic space to search for the enhanced piezoelec- 
tricity. 

The optimization problem can be summarised as to 
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TABLE I: Values of the homogenized piezoelectric strain co- 
efficients dj,/(in pC/N), piezoelectric stress coefficients Cju (in 
C/vc?) and free dielectric permittivity (in kq) of single 
crystal BaTiOg. (kq is the permittivity of free space). 



Method 


di5 


^31 


^33 


ei5 


631 


633 




T 
'=33 


Simulation 
Experiment 


560.7 
564.0 


-33.7 
-33.4 


94.0 
90.0 


34.2 
34.2 


-0.7 
-0.7 


6.7 
6.7 


4366 
4380 


132 
129 



'Ref. |3: 



find (7) that maximize, 

/(7) = ^33 

subject to, 

: — TT < 7 = {(j>,6,ip) < TT for single crystal 

: < 7 = A^e, ^^0) < 7^/2 1 nolvcrvstal 
:O<7^(fT^,^e,a.0)<5 / polycrystal J 

(8) 

An interface between our modified simulated annealing 
algorithm (implemented in Matlab) and homogenization 
program (implemented in Fortran) is created to call ^33 
(the objective function) at each iteration of the optimiza- 
tion. The temperature T is set to fall geometrically by 
20% from each of the previous step fc, i.e., Tk+i = O.STfc. 
Ideally, we must start the iteration with an initial guess of 
the design variables randomly picked up from ct G [0,5] 
and /i e [0,7r/2]. To verify the correctness of the al- 
gorithm, the optimization procedure is first applied to 
the case of single crystal BaTi03. Thus we started with 
discrete Euler angles 0, 0, and "0 alone without going 
to the assumption of distribution of grain orientations 
since a single crystal is devoid of any grain structure^^. 
All the three angles are allowed values between limits 
— TT < (0, 6, ip) < +TT. This is to accommodate the analy- 
sis of orientational dependence of piezoelectricity in sin- 
gle crystal BaTi03 as well. 

III. RESULTS AND DISCUSSION 

A. Objective: single crystal piezoelectricity 

Since the ^33 component (in tensor form 1^33 is ex- 
pressed as £33 — dsssE^, where £33 is the piezoelectric 
strain along the electric field £^3) of piezoelectric coef- 
ficient tensor better describes the piezoelectricity along 
the spontaneous polarization direction of the BaTi03'^'^, 
we designate it as the objective function for our first op- 
timization problem. First, we would check the robust- 
ness of the homogenization implementation before start- 
ing with the optimization. Good agreement between our 
homogenization results on single crystalline BaTi03 and 
the experiment by Zgonik et. ali^ is as shown in Ta- 
ble m The comparison on Table U would provide only 
a check to the correctness of the algorithm. In fact, in 




Iterations 

FIG. 3: Variation of effective piezoelectric stress coefficients 
(a) 633, (b) 631 and (c) 615 along with piezoelectric strain co- 
efficients (d) di5 and (e) dai in single crystal BaTi03 against 
iteration number. 



our simulation we use the elastic stiffness C^i, , piezoelec- 
tric stress constants e^^, and clamped dielectric permit- 
tivity nfj from the measurements of Zgonik et ali^. The 
present simulation is limited in its detail in the sense that 
we have not accounted the domain structure of the grains 
since they are assumed to be composed of as single do- 
main. Nonetheless, the extrinsic contributions associated 
with the displacement of domain walls at external fields 
has a profound influence on the dielectric, mechanical and 
piezoelectric properties of ferroelectric materials^'^'^'^. 
Moreover the size effect of grains and domains are critical 
to the piezoelectricity of ferroelectrics as well^i^. 

The evolution of the objective function c?33 with 
the temperature in single crystal BaTi03 is shown in 
Fig. [T] The piezoelectric coefficient c?33 obtained after 
optimization (^33 — 223.7 pC/N) compares well with 
experimental^ (also see. Table |TT| and other theoretical 
resultsS^. The solution 0, 9 and ip are -2.182, 0.873 and 
-0.175 radians respectively. This if expressed in degrees 
(0 = -125°, e = 50° and tP = -10°) would correspond to 
one of the < 111 > directions of the BaTi03 single crys- 
tal along which the maximum piezoelectric coefficient of 
^33 = 203 pC/N is measured2^. 
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TABLE II: Optimized piezoelectric coefficient daa along with 
experimental values of single crystal and polycrystal BaTiOa . 

Optimum objective function, d-^^ (pC/N) 
Authors Single crystal Polycrystal 

Present 223.7 " 270.7 ' 

Wada et al. " 203.0 

Zgonik et al. '' 90.0 

Bechmann " - 191.0 



"Optimization is achieved when the crystal reaches the orientation 
Euler angles (</> = -2.182, 9 = 0.873, ip = -0.175). This is a {111} 
orientation. 

'Optimal microstructure is characterised by the orientation dis- 
tribution parameters (fig = 0.785, erg = 0.1, fi^ = 1.134, cr^ = 3.3, 
p.^ = 0.873, 17^=4:) 
■^Experiment in Ref. |33| for [lll]-poled single crystal BaTiOa 
■^Experiment in Ref. 32 for [001]-poled single crystal BaTiOa 
■^Experiment in Ref. |58j for poled polycrystalline BaTiOa 

The plots on Figs. HJa) and (b) reveals that the piezo- 
electric coefficient ^33 always confined vifithin ±223.7 
pC/N. Here FZ-projection (circles) on Figs. [2] (a) and 
(b) shows peaks at 6* = ±50° and at ±140°. That 43 
here follows the same pattern shown in earlier studies^i^ 
about the rotation (9) of the single crystal provides an ad- 
ditional proof of the efficiency of the algorithm and show 
how comprehensive is the solution space. Also, it can 
be seen that the [XZ-projection (rectangles) on Figs. 
(a) and (b)] angles (j) and have no visible influence on 
the piezoelectricity of BaTi03. The (^33 exhibits both 
highs and lows irrespective of the values of (j) and ip as 
dictated by the symmetric reduction of the transforma- 
tion equation in Eq. ([1]). Also, since the spontaneous 
polarization in tetragonal BaTi03 is along the [001] di- 
rection, any rotation of that axis of the crystal affects the 
piezoelectricity of the crystal measured in a local frame 
of reference. This is evident considering the fact that 
the Euler angle 9 measures the amount of rotation of the 
[001] axis of BaTi03. (This aspect would be discussed in 
detail in the coming section). In other words, for single 
crystalline BaTi03 to display enhanced piezoelectricity, 
one should cut the single crystal at an angle away from 
the polar axis. 

The objective function ^33 attains the optimum value 
of 223.7 pC/N at a temperature of 6.4 [Fig. [T] and 
at the thirteenth iteration. (We have listed the c?33 
parameter obtained at the optimization in Table |IT] 
along with the experimental values). After the symmet- 
ric reduction of the coordinate transformation equation 
dkij = O'kpO-im.O'jndpmm whcre ttij are the elements of Eu- 
ler transformation matrix in Eq. ([T]), yields 

dgg = {d\^ + dgi) ^in^ 6* COS 6* + 43 cos^ 9 (9) 

for tetragonal ^mm symmetry of BaTi03. Here the 
primed coefficients corresponds to that of the sponta- 
neously polarized BaTi03. Almost 89 % of the contri- 
bution to the enhanced c?33 is from the first term of this 



Temperature 




Iterations 



FIG. 4: Objective function (piezoelectric coefficient ^33) at 
each iteration in FE polycrystalline BaTi03. The circles in- 
dicate the value of daa at each iteration. The line shows the 
values retained by the optimization algorithm at each temper- 
ature step. The dotted lines indicates the experimental values 
(= 203 pC/N) obtained for [lll]-oriented engineered-domain 
BaTi03 (Ref. 133^ and (= 191 pC/N) for poled polycrystalline 
BaTi03 (Ref. [53). (Here the temperature is a decreasing con- 
trol parameter of the optimization). 

equation containing the shear constant 1^15. i.e., while 
evaluating the terms in the above Eq. [SI we obtain the 
contribution from the first term as 198.8 pC/N and the 
rest of 223.7 pC/N is derived from the last term 1^33 cos"^ 9. 
We have used the simulation results of d'^^ from Table H] 
for this calculation. This corroborates the notion of the 
relation between polarization rotation, shear constant c?i5 
and the eventual piezoelectric enhancement^. 
The tetragonal phase of BaTi03 has^S 

d-i?, = 2e3isi3 + 633333 (10) 

Table Hill provides the set of piezoelectric coefficients djv, 
ejn and compliance s^^ obtained at the point where the 
1^33 attains the optimal value in BaTi03 single crystals. 
Obviously, the major contribution to the c?33 comes from 
the second product of the above equation (Eq. [TU)) . As 
one can see from Tables [T] and IIIIl the piezoelectric co- 
efficient 633 jumps from 6.7 to 27.3 C/m^. This big in- 
crease contributes the major part of the enhancement of 
^33. The variation of Cjv and dj^ are shown in Fig. [31 
Fig. [3l facilitates the understanding of the geometrical 
relations contributing to the enhancements of piezoelec- 
tricity in BaTi03. In addition to this as shown in Eq. [51 
^33 becomes a linear combination of dj^ after coordinate 
transformation. This aspect would also contributes to 
the enhancement of ^33. 

B. Objective: polycrystal piezoelectricity 

The optimization of polycrystal BaTi03 is treated in 
this section. The novelty in the present work is that 
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TABLE III: Piezoelectric coefficients dj^ (in PC/N), ejv (in C/m^) and compliance s^^ (in 10 m^/N) obtained at the point 
of optimal dsa in BaTiOa. 



Phase 


633 


631 


615 




rfl5 


Sll 


Sl2 


Sl3 


S33 


S44 


Single crystal 


27.3 


-17.7 


4.7 


-144.2 


115.4 


6.7 


-2.5 


-1.4 


5.9 


17.9 


Polycrystal 


30.0 


-13.8 


14.4 


-124.0 


254.5 


6.3 


-1.9 


-2.2 


6.9 


18.0 




FIG. 5: Variation of the objective function, the piezoelec- 
tric coefficient d33 with the design variables, mean and stan- 
dard deviation of each of the angle distribution in polycrystal 
BaTiOs. Each of the subplots shows the dss against the de- 
sign variables, (a) fie, (b) ere, (c) /x^, (d) a^, (e) fi^ and (f) 
corresponding to each iteration. 



the optimization procedure being applied to the ferro- 
electrics. Precisely, a stochastic optimization is devel- 
oped used to find out the best grain configuration which 
will output a piezoelectric polycrystal wherein the piezo- 
electricity is maximum compared to the other possible 
configurations. We analyze the most general case with 
((T0, 0-0, cr^) e [0,5] and {fi^, fig, fi^) € [0,7r/2] as given 
by Eq. [S] The results are shown in Fig. S) The initial 
and final temperature of the optimization is set to be 1 
and 1. 0634 X 10"'* respectively. Each step (there would be 
42 steps to fall into the final temperature) constitutes 24 
iterations and as a whole the final solution is realized in 
1008 iterations. The objective function converges with a 
value dss — 270.7 pC/N which is much higher than both 
[001] poled and [111] poled single crystal values (see Ta- 
ble |lTl. Here, the ^33 is obviously enhanced by a factor 
of 3 from the [001] poled single crystal value of ^33 « 90 
pC/N as is seen in Fig.|4l (Here the ^33 of textured poly- 
crystal is being compared to that of the [001] poled single 
crystal). The optimal value of d^s is also higher than the 
corresponding value of an optimally oriented single crys- 
tal (where the value was (^33 = 223.7pC/N as shown in 
Table HB . Also this is much higher than the poled poly- 
crystal experimental value of ^33 = 191 pC/N obtained 
by Bechmann^. 



o 

Q. 




FIG. 6: Variation of the homogenized piezoelectric coefficient 
^33 with the orientation distribution parameters ag and /ig 
(in degrees) for the polycrystalline BaTi03. The other two 
Euler angles {<j),^) are kept at zero.. 



The solution obtained is = 0.785, ag = 0.1, fi^ = 
1.134, (70 — 3.3, fi^ = 0.873, — 4. It means 6s (which 
measures the canting of the c-axis of the crystallographic 
grain) is kept at a small standard deviation of 0.1 but 
around a mean value of 0.785 radians (« 45°). Neverthe- 
less, the other two angles are distributed with larger stan- 
dard deviations close to the random. Thus the present 
solution suggests one should keep the Euler angles <j> and 
tp related to the orientation of ab-plane of the crystal- 
lites to be in random while the orientation 6, of c-axes 
is kept close to 45° but with a marginal standard devia- 
tion. A similar kind of result in polycrystalline BaTiOs 
was obtained numerically by Garcia et al.^P. They have 
shown that maximum piezoelectric response (c?33 or ^31) 
is exhibited by polycrystal ferroelectrics possessing spe- 
cific crystallographic textures. 

The dependance of the spread and shape of the ori- 
entation distribution on the piezoelectric properties are 
displayed in the three-dimensional plots on Figs. [S) It 
can be seen that the effective piezoelectric constant c?33 
clearly shows a distinctive dependence on the distribu- 
tion parameters pertaining to the Euler angle 6 [as can 
be seen from the pattern shown in Figs. [5] (a) and (b)]. 
(The contour on the YZ- projection on the plots shows 
the values of ^33 for the parameters fig, ag, /i^, cr^, /i^, 
cr^). The piezoelectric coefficient follows a specific pat- 
tern for both fig and ag shown on subplots Figs. [5] (a) 
and (b) respectively. ^33 shows peaks around the mean 
fig « 0.7 — 0.8 radians which would be around sa 45°. 
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Also, dss steadily increases its value as the standard de- 
viation CTg — >■ 0, which is clearly a tendency towards an 
aligned and textured ceramic material. 

Next, we will study the influence of randomness on the 
enhancement of piezoelectricity in ceramic BaTiOa. An 
important aspect observed in the optimization study is 
that in polycrystals, unlike in single crystals, the orien- 
tation of the ab-plane of the crystallites do play a role 
in determining its piezoelectricity. Another simulation is 
done to verify this point as shown in Fig. [51 Here the sim- 
ulation is done keeping the (j) and at zero while letting 
the fig and ag varying. In contrast to the results obtained 
for polycrystal without any constraint (see Figs. 2] and 
[5]), this simulation results in ^33 = 236 pC/N which is 
close to that of [lll]-oriented single crystal. This would 
further point out that c-axes of the crystallites should 
be constrained to have a specific orientation while the 
ab-plane need not be kept at a specific texture but at 
random. This condition will deliver a ceramic piezo- 
electric material possessing better piezoelectricity than 
any other phase (optimally oriented single crystal or ce- 
ramic). Hence randomness in the orientation of grains, 
if utilized judiciously, could be useful for manufacturing 
piezoelectric ceramics which outperform single crystals. 
Another point to noted here is the role played by the 
shear constant di^ in the enhancement of ^33. Here the 
analysis of the the shear constant in the enhancement 
of piezoelectric response is complicated because of the 
presence of grain boundaries in polycrystal FEs. Yet, it 
would be important to know the value of di5 to shed light 
on the origin of piezoelectric enhancement. The value of 
the shear constant obtained at the maximum of ^33 is 
di5= 254.5 pC/N (see Table As in single crystals 

the enhancement in rf33 in polycrystals is accompanied 
by a sizeable di5 value in BaTi03. A similar conclusion 
was drawn by Garcia et alM^, where they show that in 
partially textured ceramics with a di5 larger than 1^33 is 
a good choice for enhancement of polycrystal ^33. 

The optimal orientation distribution across the six 
faces of the microstructure (representative volume ele- 
ment) is shown in Fig. [T] Since the angle 9 influences 
the d^s the most, that only is shown. It is seen that 
(Fig. [21) the angles are mostly clustered around 0.8 ra- 
dians (~ 45°). The effective piezoelectric coefficient is 
dictated by the transformation of coordinates involving 
the Euler angles [see Eqs. ([T]) and ([5])] and hence is in- 
fluenced by the rigorous combinations provided by the 
randomness. One of the possible reasons for higher piezo- 
electricity in FE ceramics with a certain pattern of grain 
distribution is the role played by the intrinsic polariza- 
tion of the grains. The misalignment of the polarization 
in the neighboring grains would impart residual stresses, 
which couples with the electrical field to enhance the po- 
larization and thereby the piezoelectricity. In the optimal 
distribution, the polarizations of the juxtaposed grains 
are configured in such a way that they will contribute 
positively to the overall increase of piezoelectricity of the 
polycrystal. 




FIG. 7: The distribution of Euler angle 9 in optimized poly- 
crystalline BaTiOs. Here the angles in radians are mapped 
in colors. The six planes shown from (a)-(f) are (100), (001), 
(Too), (010), (OTO) and (GOT) respectively. 



In summary, we have introduced a global optimiza- 
tion technique to identify the ideal configuration of both 
single and polycrystal FE. This method can be used to 
the FEs of any symmetry and to any objective func- 
tion involving electromechanical coefficients. Here this 
method is applied to the perovskite, tetragonal ferroelec- 
tric BaTi03. The paper predicts configuration of crys- 
tallites which maximizes the piezoelectricity of the poly- 
crystalline ferroelectrics. The orientation distribution re- 
sulted from electrical and mechanical poling is limited 
by many factor The alignment of domains with the 
poling field is by switching its polarization to an equi- 
librium position defined by the symmetry. The as-grown 
polycrystal before poling has domains randomly oriented. 
Thus the poling doesn't yield the polarization vectors 
(of domains) all aligned along the field since there may 
not exist an equilibrium orientation for some domains 
along the poling field direction. Apart from being pro- 
vide a guide to experiment, essentially it doesn't suggest 
a particular method (for instance, mechanical or elec- 
trical poling, tape-casting, templated grain growth, hot 
forging etc.) to realise the predicted crystallite configu- 
ration. A plausible approach is to adopt a combination of 
suitable processing method (to achieve a certain degree 
of preferred orientation in the unpoled state) and poling. 

In this work, we have observed significantly enhanced 
piezoelectric response in ferroelectric polycrystals at cer- 
tain grain configurations. In single crystal BaTiOa the 
piezoelectricity is found to be larger along a nonpolar di- 
rection away from the polar axis. We have optimized the 
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ferroelectric ceramic by design at the microstructure level 
for piezoelectric applications. The solution obtained from 
the optimization procedure results in a three-fold en- 
hancement of piezoelectricity in the ceramic phase com- 
pared to the single crystalline phase. If we use the ran- 
domness of the grain orientations judiciously the ceramic 
can replace even the oriented single crystals in piezoelec- 
tricity. A plausible reason behind the anisotropy shown 
by both rotated single crystal and polycrystal FE could 
be the macroscopic symmetry. The crystallographic sym- 
metry is characterised by the anisotropic (or isotropic) 
physical properties. The product phase obtained after 
the crystal orientation might be different from the par- 
ent phase in symmetry. The insight obtained from the 



optimization have the potential to advance the design 
and discovery of complex FE configurations with supe- 
rior piezoelectric performance. Further studies in this 
direction in relaxor ferroelectrics, where the single crys- 
tals display larger piezoelectricity, could inaugurate new 
possibilities in technological applications which involves 
the requirement of high piezoelectricity. 
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